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RECENT BES RESULTS ON ^(3770) AND D MESON PRODUCTION AND 

DECAY* 

G. RONG 

(for the BES Collaboration) 
Institute of High Energy Physics, Beijing, 100039, China 

Using a data sample of 17.7 pb" 1 collected at 3.773 GeV with the BES-II detector at the BEPC, 
the cross sections for D°D and D + D~ productions at 3.773 GeV have been measured. From 
the data sample about 33 pb" 1 taken around 3.773 GeV, 7696 ± 199 ± 369 and 5381 ± 128 ± 188 
singly-tagged neutral and charged D mesons are accumulated, respectively. In the system recoiling 
against the singly-tagged charged D sample, 3 purely leptonic decay events of D + — ► jj, + v are 
observed, which yields the branching fraction to be BF(D + -> /j+v) = (0.120+oog3+°oJg)% and 
the decay constant /d = (365l 113 l 2 8) MeV. From the singly tagged D sample, the semileptonic 
decay branching fractions for D° — » K~e + v, D° — > ir~e + v and D + — + K e + v are measured to 
be BF(D° -> K~e + u e ) = (3.52 ± 0.36 ± 0.25)%, BF(D° -> 7r"e + !/ e ) = (0.36 ± 0.14 ± 0.03)% 
and BF(D + -> K°e + u e ) = (8.64 ± 1.51 ± 0.72)%. The vector form factors are determined to be 
|/+ (0)| = 0.75 ± 0.04 ± 0.03 and |/+(0)j = 0.76 ± 0.15 ± 0.06. The ratio of the two form factors 
is extracted to be |/£(0)//+ (0)| = 1.01 ± 0.20 ± 0.08. The ratio of the decay widths is measured 

to be r(D° -> K~e + v)/Y(D + -» K°e + v) = 1.04 ±_0.21 ± 0.08. From the data sample of about 27 
pb -1 , the evidence of 1^(3770) — » J/ipir + -K~ non-DD decay is observed. The branching fraction is 
determined to be 5^(^(3770) -» J/?/>7r + 7r") = (0.34 ± 0.14 ± 0.08)%, corresponding to the partial 
width of T(tp(3770) -* J/^tt+tt") = (80 ± 32 ± 21) keV. 



I. INTRODUCTION 



At the center-of-mass energy of 3.773 GeV, the ?/>(3770) 
resonance is produced in e + e" annihilation, and the open 
charm pairs of D°D° and D + D~ are mainly produced 
from ip(3770) decays. Taking the advantage of the DD 
production, we can use the single tag method to measure 
the cross sections for D°D , D + D~ and DD productions 
at the energy of 3.773 GeV; using the double tag method 
we can measure some absolute decay branching fractions 
of neutral and charged D mesons. 



Width nD°^K-e +v) 
r(D+— K e+v) 

to a i\oi\-DD final state 



and the evidence of ^(3770) decay 



In this paper, we report some preliminary results on 
measurement of the cross section for DD production, 
measurement of the decay constant fn, measurements 
of the form factors /+(0) and f+(0) in the semileptonic 
decays of D° meson, measurement of the ratio of partial 



?S II. CROSS SECTION FOR DD PRODUCTION 



The data used in this analysis were collected at the 
center-of-mass energy of 3.773 GeV with the Beijing 
Spectrometer p| at the Beijing Electron Positron Col- 
lider. The total integrated luminosity of the data is 17.7 
pb" 1 . The measurements of the cross sections for the 
D°D°, D + D~ and DD productions are made based on 
the analysis of singly tagged D° and D + samples. At 
center-of-mass energy y/s — 3.773 GeV, the D° (Through 
this paper, charge conjugation is implied.) and D + arc 
produced in pair via the process of 



e ' e 



D°D°,D + D' 



(1) 



*The talk given at 39th Recontres de Moriond on Electroweak In- 
teraction and Unified Theories, LaThuile, Aosta Valley, Italy, 21-28 
March, 2004. 



The totally observed number N D o a {N D + ) of D° {D+) 
mesons and the observed cross sections <7° b ,f— (ct^+ d _) 



D°D 



are related by a relation 



D°D° and DD productions to be 



obs 

d°d" 



and 



N. 



LK' 



2 x L x B x e 



(2) 



cr£%- = (2.52 ±0.07 ±0.24) nb, 



ct d°d° = ( 3 - 26 ± °- 09 ± °- 25 ) nb - 



„obs 
a D+D- 



N 



D 



2 x L x B x e 



(3) 



where L is the integrated luminosity of the data set 
used in the analysis, B is the well-known branching frac- 
tions 0] for decay modes in question, and e is the Monte 
Carlo efficiency for reconstruction of this decay mode in 
the data analysis. 

In the measurements of the cross sections, the singly 
tagged neutral and charged D meson are made by fully 

reconstructing one D meson from the D D (D + D~) 
pair in the invariant mass spectra of the daughter parti- 
cles from the D decay. Taking the advantage of the DD 
pair production, we use kinematic fit to some specific par- 
ticle combinations to improve the ratio of signal to noise 
and mass resolution in the invariant mass spectrum. The 
distributions in the fitted masses of Krmr (m = 1, or 2, 
or 3) combinations, which are calculated using the fit- 
ted momentum vectors from the kinematic fit, are shown 
in Fig. ^ The signals for neutral and charged D mesons 
production are clearly observed in the fitted mass spectra 
as shown in the Fig. ^ 
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FIG. 1: Distribution of the fitted masses of the Kmir (m=l, 
or 2, or 3) combinations for three singly tagged modes, where 
Fig. (a) and Fig. (b) are for the decay modes of D° — » K~~-k + 
and D° — ► K~tt + -k + -k~ , respectively, and the Fig. (c) is for 
the decay mode of D + — ► K~-k + tt + . 

We measure the observed cross sections for D + D~, 



a<g s B = (5.78 ± 0.11 ± 0.45) nb. 

After correcting to the effect of Initial State Radiation 
and the vacuum polarization, we obtain the tree level 
cross sections to be, 



tree 
T D+D- 



tree_ 
D°D° 



(3.23 ±0.09 ±0.34) nb, 



(4.18 ±0.12 ±0.37) nb. 



III. 



a% e £ = (7.42 ± 0.14 ± 0.67) nb. 



PURELY LEPTONIC DECAYS OF D+ 
AND DECAY CONSTANT f D 



fl u 



K u tx- 



tt + modes with sub-resonances 



The D mesons are reconstructed in K 

K°K~, K + R-TT-, K°TT--r--r + 

K + 1T + TT~1T~TT~ and7T~7r" 

decaying as Kg — ► ir + ir~, and ir° — > 7 7. Fig. shows 
the distributions in the fitted masses of mKnir (m =1,2 
and the n = 1,2,3) combinations. The signals for D~ 
production are clearly observed in the fitted mass spectra 
as shown in the Fig. We accumulate 5381 ± 128 ± 188 
singly tagged D" mesons in total. 

The system recoiling against the singly tagged D~ in 
each of the events as shown in the Fig. 2 are examined 
for consistency with the decay D + — > n + v. It is required 
that there is a singly charged track and | cos0| < 0.68, 
where the 9 is the polar angle of the charged track. There 
must be hits in the muon counter which are well associate 
within ±4a with the tracks in transverse projection and 
Z direction; the required number of hits is momentum 
dependent. 

For the candidate events, no isolated photon is allowed 
to be present, where the isolated photon is defined as 
an electromagnetic shower with the energy being greater 
than 100 MeV and the direction of the shower develop- 
ment separated by at least 20 degrees from the direction 
of the nearest charged track. Since there is a missing 
neutrino in the purely leptonic decay event, the event 
should be characteristic with missing energy E m i SS and 
missing momentum P m iss carried by the neutrino. For 
the typical purely leptonic decay events, the difference 



between the E„ 



and the P„ 



should be around zero. 



We define the deference to be 



U„ 



To select the purely leptonic decay events from the singly 
tagged D~ sample, we require the U m i SS of the candidate 
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FIG. 2: Distribution of the fitted masses of the nKmn (m=0 
or 1 or 2, n=l or 2 or 3 or 4) combinations for the 9 
single tag modes; (a), (b), Jc), (d), (e),Jf), (g), (h) and 
(i) are for the modes of D~ — > K + tt~tt~, D~ — » K°n~, 
23" -> A"°A"-, £)- -> K"+A"-7r-, L>" -> A" 7r-7r-7r + , 

and D~ — > 7r _ 7r _ 7r + respectively; (j) is the fitted masses of 
the nKrnir combinations for the 9 modes together. 



events to be in the ±3.0tJr/ mi33 region. We also require 
the momentum of the muon to be within the range of 
0.780 GeV/c to 1.115 GeV/c. Three candidate events 
satisfy the selection criteria of the purely leptonic de- 
cay events. A detailed Monte Carlo study shows that 
there are 0.31 ± 0.16 background events in the 3 can- 
didate events. The efficiency to reconstruct the purely 
leptonic decay events is 0.405 ± 0.011. These result the 
purely leptonic branching fraction to be 



BF{D^ 



/' 



(0.120; 



-0. 092+0. 010\ 
-0.063-0. 009 ) 



and decay constant of 

f D = (365±g|±^) MeV, 

where the first errors are statistical and second system- 
atic which arise from the uncertainties in the measured 
branching fraction, the CKM matrix element \V c d\ and 
the lifetime of D + : the systematic uncertainties in the 
measured branching fraction arise mainly from the un- 
certainties in the [i + identification, tracking efficiency, 
background subtraction and U m i SS cut. The central value 
of the Sd is consistent with that measured with BES-I 
detector HI. 



SEMILEPTONIC DECAYS OF CHARGED 
AND NEUTRAL D MESONS 



Branching fractions for D° 

D°->TT-e + v 



K e v and 



-r=rO 



The singly tagged D mesons are obtained by recon- 
structing the four decay modes of D — > K + n~ , D —* 



D° 



K+TT-TT-TT+, D -> K°n+n- and D" K^ 

Fig. |21 shows the fitted mass distributions of Knir (n = 
1, 2, 3) combinations. We totally accumulate 7696±199± 
369 singly tagged neutral D meson sample. Fig. 0Ja) 
and Fig. EJb) show the distribution of the fitted masses 
of the Kwk combinations for the events for which the 
D° — > K~e + v e and D° — > ir~e + v e candidate events 
are observed in the system recoiling against the singly 
tagged D . After subtracting numbers of background 
events (sideband background events and the background 
events due to other semileptonic or hadronic decays), 
100.5±10.2 and 10.3±3.9 signal events for D° -> K~e + v e 



and£> ( 



o 



decays are retained. We measure the 



branching fractions to be 

BF{D° -+ K-e+v e ) = (3.52 ± 0.36 ± 0.25)%, 



BF(D° 



(0.36 ±0.14 ±0.03)%. 
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FIG. 3: Distributions of the fitted invariant masses of (a) 
K+tv-, (b) K+-K-7V+TT-, (c) K^ir + ir~ and (d) K+tx'-k com- 
binations. 
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FIG. 4: Distributions of the fitted invariant masses of Knn 
combinations for the events for which (a) the D° — > K~e + u e 
and (b) the D° — > 7r~e + z/ e candidate events are observed in 
the system recoiling against the tagged D . 



B. Branching fraction for D~* 



K e 



Candidate events D H 



*V 



is selected from 
the surviving tracks in the system recoiling against the 

tagged D~ . To select the D + — > X e + i/ e events, it is 
required that there are only three charged tracks, one of 
which is identified as an electron and the other two tracks 
as 7r + and 7r~ . The candidate events are required to sat- 
isfy the requirement |Z7 TO j ss | < 3<7;/ misl> , where the (Tu mUa 
is the standard deviation of the U m iss distribution. 

Fig. [5] show the distribution of the fitted invariant 
masses of the Knit combinations for the events for which 
the D + — > K e + v e candidate events are observed in the 
system recoiling against the singly tagged D~ . We mea- 
sure the branching fraction to be 



BF{D^ 



ifV 



5 ) = (8.64 ±1.51 ±0.72)%. 
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FIG. 5: Distributions of the fitted masses of Knit combina- 
tions for the events for which the D + — » K e + v e candidates 
are observed in the system recoiling against the tagged D~ . 



The form factors |/^(0)| and |/+(0)| can be extracted 
by using the measured values of the branching fractions 
and the lifetime of the D° meson. Inserting the values 
of |V C8 | = 0.996 ± 0.013, \V cd \ = 0.224 ± 0.016 and the 
lifetime r", = (411.7 ± 2.7) x 10" 15 @ into equation © 
and |JSJ, the form factors are extracted to be |/^(0)| = 
0.75 ±0.04 ±0.03, |/^(0)| = 0.76 ±0.15 ±0.06, where the 
first errors are statistical and the second are systematic. 
We extract the ratio of the two form factors to be 



|/+(0)//f (0)| = 1.01 ± 0.20 ± 0.08. 



D. 



The ratio of T(D° -> K-e+u)/T(D + -> K e+u) 



Isospin symmetry predicts that the ratio of partial 
widths T{D° -> K-e+v)/T(D+ -> K°e+v) should be 
unity. The ratio can be obtained from the measured 
branching fractions for D° — > K~e + u and D + — > K e + v, 
and the lifetimes of the D° and D + mesons. 



r(L>° -^ K-e+v) _ BF(D° -> K~e + v) t d + 
T(D+ -> K°e+v) BF(D+ -» ife+v) t d« 



(6) 



C. Form factors /£(0) and /f (0) 

The relations Q| [5| between the decay widths and the 
form factors are, 

T(D° - K-e+v e ) = 1.53 |^ cs | 2 |/f (0)| 2 x lO 11 ^ 1 , (4) 



T(D° -> n-e + v e ) = 3.01 |K d | 2 |/;(0)| 2 x lO 11 ^ 1 . (5) 



where the t d + and t d o are the lifetimes of D + and D° 
mesons. Inserting the measured branching fractions, and 
the lifetimes of the D + and D° mesons, we obtained the 
ratio to be 



T(D° -» K-e+v) 
Y(D+ -* K°e+p) 



1.04 ±0.21 ±0.08, 



which is consistent with that the isospin conservation 
held in D meson decays. 



V. EVIDENCE OF V(3770) -> J/tpir + ir~ 

The 0(3770) resonance is believed to be a mixture of 
the l 3 Di and 2 3 S' 1 states of the cc system (g. Since its 
mass is above open charm-pair threshold and its width is 
two orders of magnitude larger than that of the 4>(2S), it 
is thought to decay almost entirely to pure DD .7] . How- 
ever, recently some theoretical calculations point out that 
the ■0(3770) could decay to non— DD final states [8| Q- 

To search for the decay of 0(3770) — > J/i(jn + n~ , 
J/ip — > e + e~ or /i + /i~, /i + /i~7r + 7r~ and e + e~7r + 7r _ 
candidate events are selected. These are required to 
have four charged tracks with zero total charge. Can- 
didate events are subjected to four-constraint kinematic 
fits to either the e + e~ — ■> /i + fi~ ir + ir~ or the e + e~ — > 
e + e~7r + 7r~ hypothesis. 

Fig. ICTa) shows the dilepton masses determined from 
the fitted lepton momenta of the accepted events. There 
are clearly two peaks. The lower mass peak is mostly due 
to 0(3770) — ► J/-07r + 7r~, while the higher one is due to 
ip(2S) — > J/-07T+7T - . Since the higher peak is produced 
by radiative return to the tp(2S) peak, its energy will be 
approximately 3.686 GeV, while the cm. energy is set to 
the nominal energy in the kinematic fitting. Therefore, 
the dilepton masses calculated based on the fitted lepton 



momenta from ij){2S) — > J/tJjtt + tt , J/ip 



l+l- 



shifted upward to about 3.18 GeV. The fit to this peak 
yields a J/ip signal of 17.8 ± 4.8 events. 

There is a contribution from ip(2S) — > J/tJjtt + tt~ pro- 
duced by radiative return to the tail of the 0(25) that can 
pass the event selection criteria and yield fitted dilepton 
masses around 3.097 GeV. This is the main background 
to '0(3770) — y J/07r + 7r~, as shown in Fig. [B£b). Here 
the histogram shows the dilepton mass distribution for 
-0(25*) — > J/tpir + ir~~ from our Monte Carlo simulation. 
The higher peak is due to the radiative return to the 
ip{2S) peak, and the lower peak is from radiative return 
to the tail of the ip(2S). The points with error bars show 
the total contribution from ip(2S) and 0(3770) produc- 
tion and decay. From the simulation, we estimate that 

6.0 ± 0.5 ± 0.6 out of 17.8 ± 4.8 events in the peak near 

3.1 GeV in Fig. EJa) are due to ip{2S) -> J/07T+7T-, 
where the first error is statistical and the second is the 
systematic arising from the uncertainty in the ip(2S) res- 
onance parameters. The probability that the 17.8 events 
observed are due to a fluctuation of the 6.0 events is 
3.8 x 10" 4 . 

After background subtraction, 11.8 ±4.8 signal events 
are remained. The branching fraction for the non— DD 
decay 0(3770) — > J/-07r + 7r~ is measured to be 

BF(t/>(3770) -> J/Vtt+tO = (0.338 ± 0.137 ± 0.082)%. 



Using r to t from the PDG Q{, this branching fraction cor- 
responds to a partial width of 



T(^(3770) -> J/iP-k+tt-) = (80 ± 32 ± 21) keV. 



More detail may be found in Ref. [T3| and Ref. [Tlj . 
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FIG. 6: Distributions of dilepton masses for (a) data and 
(b) Monte Carlo sample. The hatched histogram in (a) is 
for J/ip — » fj, + n~ , while the open one is for J/ip — » e + e~. 
The histogram in (b) is for ip(2S) — » J/ipir + ir~ , while the 
points with error bars are the sum of 0(3770) — > J/ipiv + iT~ 
and 0(2S) — ► J/ipir + ir~ . (c) Distribution of mass recoiling 
against the ir + n~ system calculated using measured momenta 
for events that pass the kinematic fit requirement, where the 
hatched histogram is for J/ip — » fi + fi~ and the open one is 
for J/ip — + e + e~ . 
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